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ABSTRACT 

We present a detailed study of an estimator of the Hi column density, based on a combina- 
tion of Hi 2 1cm absorption and Hi 2 1cm emission spectroscopy. This "isothermal" estimate is 
given by A^hi.iso = 1-823 x 10'** / [rtot x Tb] / [1 - e"^""] dV, where Tto, is the total Hi21cm 
optical depth along the sightline and Tb is the measured brightness temperature. We have used 
a Monte Carlo simulation to quantify the accuracy of the isothermal estimate by comparing 
the derived A^hi.iso with the true Hi column density A^hi- The simulation was carried out for a 
wide range of sightlines, including gas in different temperature phases and random locations 
along the path. We find that the results are statistically insensitive to the assumed gas tem- 
perature distribution and the positions of different phases along the line of sight. The median 
value of the ratio of the true Hi column density to the isothermal estimate, A^hi/A^hi,iso, is 
within a factor of 2 of unity while the 68.2% confidence intervals are within a factor of 3 of 
unity, out to high Hi column densities, < 5 x lO^-* cm"^ per 1 km s ' channel, and high total 
optical depths, < 1000. The isothermal estimator thus provides a significantly better measure 
of the Hi column density than other methods, within a factor of a few of the true value even 
at the highest columns, and should allow us to directly probe the existence of high Hi column 
density gas in the Milky Way. 
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1 INTRODUCTION 

The neutral atomic hydrogen column density A'hi is an important 
input to our understanding of the interstellar medium (ISM). For 
example, it determines whether the gas is predominantly ionized 
(for A^Hi 5 10" cm"^, in the intergalactic medium), predominantly 
neutral (for A'hi <; 2 x 10^" cm"^, in typical gas clouds in galax- 
ies), or predominantly molecular (for A'hi <: 10^- cm"^, in compact 
molecular clouds). It serves as the reference for estimates of vari- 
ous interesting quantities such as gas metallicities and abundances, 
is required to derive the gas spin temperature, and is the basic input 
for models of gas clouds. Accurate estimates of A'hi are thus critical 
for ISM studies. 

There are two standard approaches towards measuring A'hi 
in Galactic clouds. The first is based on absorption spectroscopy 
of stars and quasars in the Lyman-o- line, which develops wide 
Lorentzian wings for typical sightlines through the Milky Way (or 
external galaxies), and whose equivalent width is directly related to 
the Hi column density. Such damping wings are easily detectable 
with modem optical spectrographs for A'hi <: 10" cm"^, and offer 
accurate A'hi measurements for sightlines that do not contain much 
dust. Unfortunately, the presence of significant amounts of dust 
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along a sightline causes obscuration of the background star/quasar 
and, further, the amount of dust obscuration correlates with the to- 
tal hydrogen column density in the Milky Way. As a result, it is 
very difficult to use Lyman-a spectroscopy to measure A'hi along 
high column density Galactic sightlines, with A'hi s> 10^' cm"^. 

The second approach to A'hi measurements is via Hi 21cm 
emission studies, which directly measure the brightness tempera- 
ture Tb of the emission. For optically-thin Hi 21cm emission, the 
Hi column density is proportional to Tb, even when the emis- 
sion arises from multiple gas "clouds" with different temperatures. 
The advantage of this method is that it is easy to detect Hi 21cm 
emission along any Galactic sightline with today's telescopes. Fur- 
ther, significant progress has recently been made in correcting 
for stray radiation, received thro ugh the telescope sidelobes (e.g. 
iKalberla et"aLll2005l : ISaiaia et aU.2005 ). At low to moderate col- 
umn densities, A^hi < 10^' cm"^, comparisons between A'hi es- 
timates from the Lyman-a absorption and Hi 21cm emission ap- 
proa ches have typically yielded excellent agreemen t, to within 10% 
(e.g. lPickev & Lockmar]|l990l : IWakker et al.ll201lh . 

Unfortunately, the relation between Tb and A^m is not straight- 
forward for the general case of arbitrary Hi 21cm optical depth. 
When the optical depth is significant, one has to know both the 
location of different emitting components along the sightline and 
their individual optical depths and spin temperatures to infer A'hi 
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from the measured Tb. Assuming that the gas is optically thin only 
yields a lower limit on the Hi column density. 

Sightlines with high Hi column densities are also the ones 
that tend to have high Hi 21cm optical depths. It is thus diffi- 
cult to accurately estimate A^hi for such sightlines using either 
Lyman-a absorption or Hi 21cm emission spectroscopy. While the 
maximum Hi column density obtained in the Leiden-Argentine- 
Bonn (LAB) survey (as suming optically-thin Hi21cm emission) 
was A^Hi ~ 2x 10^^ cm"^ dKalberla et alj2005t) . significantly higher 
A^Hi values have been inferred fro m recent modelling of Hi 21cm 
emission data of external galaxies feraun et al. mPEni^ 
as well as Lyman-or absorption studies of hig h-redshift (and low- 
metallicity) gamma ray bursts and quasars dFvnbo et al.l l2009l : 
iNoterdaeme et al.ll2012h . This raises the question of whether the 
gas is indee d predominan tly molecular at A'hi > 10^^ cm"^ in the 
Galaxy (e.g. lSchavell200lh or whether atomic hydrogen can exist at 
significantly higher Hi column densities. In this Letter, we propose 
a different approach to determine both the Hi column density along 
a sightline and the error on the measurement, based on a combina- 
tion of Hi 2 1 cm emission and Hi 2 1 cm absorption spectroscopy. 



2 THE FORMALISM 

For the Hi 21cm line, the two observables are the Hi 21cm bright- 
ness temperature Tb, measured from emission spectroscopy, and 
the Hi 21cm optical depth t, measured from absorption studies to- 
wards background radio continuum sources. For a single homoge- 
nous Hi cloud, the observed brightness temperature is given by[[| 



Tb =T, x[l-exp(-T)] 



(1) 



while the Hi column density A'hi, the Hi 21cm optical depth t and 
the spin temperature Tj are related by the expression 



1.823 X 10'** X 



TdV, 



(2) 



where A'hi is in cm"-, the spin and brightness temperatures are in K, 
and dV is in km s"' , with the integral over the line profile. Note that 
the above equations make no approximations, except that T^ >> 
hv2icm/kB ~ 0.07 K (e.g. Field 1958), which should be valid in all 
astrophysical circumstances. 

For multiple Hi clouds along a sightline, equation (|2]l remains 
unchanged, except that Ts is then the column-density-weighted har- 
monic mean of the spin temperatures of the different clouds along 
the sightline; we will denote this harmonic mean spin temperature 
as (TJ. However, th e expression for Tb is mu ch more complicated 
in this situation (e.g. iHeiles & Trolan j|2003ah : 

Tb = Sf/o' Ts, [1 - exp(-r,)] x exp([-S^'^-'T,]) , (3) 

for N "clouds" with different spin temperatures and optical depths 
along the sightline, and M, clouds between us and the ;'th cloud. 

If the Hi 21cm absorption is optically thin (i.e. peak optical 
depth « I), the above expressions can be combined to obtain 



A'Hi.ar = 1.823 X 10" x TBdV 



(4) 



and one can estimate the Hi column density directly from the 
Hi 2 1cm emission spectrum. However, in the general case of ar- 
bitrary optical depth, it is not possible to determine A'hi, even on 
combining the Hi 2 1cm absorption and emission spectra. In such a 



Note that all quantities in this section are functions of velocity. 



situation, one would, in order to determine A'hi, need to know both 
the parameters (Tj, t) of individual clouds as well as the spatial dis- 
tribution of these clouds along the sightline; the latter is especially 
difficult to ascertain observationally. 

Our aim is to estimate the Hi column density, given measure- 
ments of both the brightness temperature and the Hi 2 1 cm optical 
depth. For this purpose, we define a quantity Ts,eff , akin to the spin 
temperature, by the relation 



Tb = Ts.eff x [I - exp(-Tto,)] 



(5) 



where Tb is the observed brightness temperature and t,o( is the 
total optical depth. We can, without loss of generality, relate the 
harmonic-mean spin temperature (TJ along a sightline to Tjeif by 



(Ts>=/(TB,T„)xT,,e| 



(6) 



where /(TB,rtot) is some unknown function of Tb and Ttot (and 
which also depends on the spatial distribution of clouds along the 
sightline). On replacing in equation (|2j for A'hi, we obtain 

Ttot X Tb 



A'hi = 1.823 x 10' 



/(TB,r,„,)x ■ 



rdV 



(7) 



[1 -exp(-rto,)] 

Knowledge of the function /(Tb, t,o,) would allow us to infer the Hi 
column density from measurements of Tb and Tiot. Of course, this 
function depends on the details of the sightline. However, based on 
the simulations of the next section, we find that the median value 
of the function is k 1, even in the extreme cases of high Ttot and 
high A'hi, with a spread of only a factor of a few around the central 
value. We hence set / » I to obtain 

xTb 



A', 



1.823 X 10' 



rdV 



(8) 



exp(-rtot)] 

The above e quation to estimate the H i column density was earlier 
proposed bv lPickev & Ben s^ ( 1982) , but has no t received much 
attention in the literature. IPic kev & BensonI ( Il982') refer to A'hi.iso 
as the "isothermal" estimate of the Hi column density, since equa- 
tion[8]is the same as the expression for the Hi column density when 
the Hi 2 1 cm absorption and emission arise in a single cloud with 
a fixed spin temperature. We will continue to use this terminology, 
but note that it can also be regarded as a "thin-screen" estimator 
as it only uses the total optical depth and is independent of the 
spatial disposition of the clouds. In effect, knowledge of the to- 
tal Hi 2 1 cm optical depth and the total brightness temperature Tb 
along the sightline allow one to estimate the Hi column density. 
Note that at large t, equation ((Sj reduces to 



A'h 



1.823 X 10"* (Tt„t xTB)dV . 



(9) 



3 MONTE CARLO SIMULATIONS 

In this section, we use a Monte Carlo procedure to validate the 
"isothermal" estimator by measuring the difference between the 
true Hi column density and the Hi column densities inferred from 
the isothermal (and optically-thin) methods. The approach taken is 
the converse of the situation in the real world, where we have the 
measured brightness temperature and Hi 21cm absorption profiles, 
and would like to infer the Hi column density. Instead, for the pur- 
poses of the simulation, we will assume that a sightline contains 
some fiducial "true" Hi column density A'hi in a narrow velocity 
channel, that produces both Hi 2 1 cm emission and absorption. We 
further assume that the Hi 2Icm emission and absorption in this ve- 
locity bin arise from the superposition of contributions from differ- 
ent temperature phases along the sightline (each, of course, with an 
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Hi column density lower than the total Hi column). We distribute 
the total Hi column randomly between the different temperature 
phases, and, for each distribution, compute the observed bright- 
ness temperature and the observed total optical depth. The inferred 
Tb and Ttot values are then used to determine the Hi column den- 
sity from the isothermal and optically-thin methods, and the results 
compared to the known Hi column density. This procedure is car- 
ried out for a large number of possible sightlines, covering a wide 
range of both total Hi column densities and distributions of the Hi 
column between different temperature phases. 

Although we are mainly interested in the high Hi column den- 
sity range, A'hi s> 10"' cm"^, the simulation has been carried out 
for total A^Hi values in the range 10"" cm"^ < A'hi < 10''* cm"'. The 
upper end of the range was chosen so as to probe the extremely 
high Hi columns (> 10-' cm"-) whose presence has been suggested 
by modell ing studies of Hi in local galaxies dSraun et al.l i2009l : 
iBraunI 12012) . The channel width is assumed to be 1 kms , sim- 
ilar to t he velocity resolution of the LAB all-sky Hi 2 1cm emission 
survey jKalberla et al.ll200^ . This narrow velocity width was cho- 
sen to match the narrowest full- width- at-half-maximum of Hi 2 1cm 
lines, for purely thermal line broadening of gas at a kinetic tem- 
perature of « 20 K, so that there is no loss of information due to 
under-sampling of the line profiles. Note that Hi column densities 
> 10^' cm"^ are extremely large for such a small velocity range; 
we will discuss the effect of this assumption later. 

We also note that the simulation explicity deals with the Hi 
column density within a single observed 1 kms"' velocity chan- 
nel, instead of over some physically motivated velocity width. It 
is well known that Hi in the ISM has velocity structure (i.e. cor- 
relations between neighbouring velocity channels), due to the gas 
temperature, turbulent motions, bulk motions, etc., and the net Hi 
column density along a sightline is inferred from the full Hi 2 1cm 
line profile, not merely from individual velocity channels. Our sim- 
ulation makes no assumptions about the velocity structure, because 
quantities like the gas temperature, turbulent motions and veloc- 
ity gradients are, in general, poorly known, and assumptions about 
these quantities could bias the results. Our approach of using nar- 
row (< 1 km s" ' ) channels does not lead to any loss of generality, as 
the channel width is narrow enough to properly sample the velocity 
profile. As discussed in more detail later, the primary consequence 
of not using the velocity structure in the line profile is that the error 
bars that we advocate are conservative; the advantage is that no as- 
sumptions are needed about quantities that are difficult to measure. 

The next step is to specify the spin temperatures of the differ- 
ent gas clouds along the sightline. In the classic models of the ISM, 
Hi is expected to stab ly exist in two ph a ses in pressure equ ilibrium 
with each other (e.g. iField et"Zlll969l ; Iwolfire et al.ll2003 ); these 
are the cold neutral medium (CNM, with 40 K < T^ < 200 K) and 
the warm neutral medium (WNM, with 5000 K < Ti, < 8000 K). 
Recently, there has been evidence that significant fractions of Hi 
may be in a thermally unstable phase, with 200 K ^ T^ < 5000 K 
(e.g. iHeiles & Trolandl l2003al ; iKanekar et all l2003h and that the 
CNM temperature range extends to » 20 K iHeiles & Trolanj 
l2003bh . We hence allow for gas in three temperature phases, 
the CNM, the WNM and the thermally unstable neutral medium 
(UNM), and will also use the temperature range 20 K < Ti; < 200 K 
for the CNM. As will be seen, the precise choice of the CNM, 
WNM and UNM temperature ranges do not significantly affect our 
results. 

To estimate the brightness temperature and the Hi 21cm op- 
tical depth, we require the spin temperature, not the kinetic tem- 
perature. For the CNM, the Hi 21cm transition is expected to be 



thermalized by collisions, causing Ts « T^; we will hence use 
the spin temperature range 20 K < T^ < 200 K for the CNM. In 
the case of the WNM (and, possibly, the UNM), the low particle 
number density means that collisions cannot thermalize the line; 
the spin temperature is expected to typically be lower than the ki- 
netic temperature in this phase unless there are sufficient Lyman-o- 
radiation in the cloud for resonant scattering o f these photons to 
coupl e the spin and kinetic temperatures (e.g. lField|[l958l ; iLisztl 
l200lh . We have hence assumed tha t T^ < Ti; in the WNM, and 
have used the results of lLiszll ( 1200 ll) to relate the spin temperature 
to the kinetic temperature; this essentially meant using the range 
2000 K < Ts < 5000 K for the WNM. The UNM spin tempera- 
ture range is set to 200 K < Tj < 2000 K, intermediate between 
the CNM and WNM phases. Note that this procedure effectively 
allows for gas at all spin temperatures 20 < Tj < 5000 K. 

In order to divide the total Hi column density into different 
components along the sightline, we first fix the fraction of gas in 
the three temperature phases. The simulations were carried out for 
the following six cases: 

(i) Half the gas is in the CNM, and half in the WNM, 

(ii) Each of the three phases contains one-third of the gas, 

(iii) 50% CNM, 40% WNM and 10% UNM, 

(iv) 50% CNM, 10% WNM and 40% UNM, 

(v) 90% CNM and 10% WNM, and 

(vi) The fraction in each phase varies randomly from run to run. 

Having fixed the total Hi column density and the gas fraction 
in each phase, we know the Hi column density in each phase. We al- 
low for the possibility that there could be multiple "clouds" from a 
given phase that contribute to the emission/absorption in the veloc- 
ity channel by allowing gas in each phase to be further sub-divided 
into multiple components. We then randomly assign a spin tem- 
perature to the first component of each phase, subject to the con- 
straint that it lies within the Ts range of the phase. The Hi column 
density of the component is also randomly chosen from the range 
I X 10" cm"- to 5 X 10^^ cm"^. We then add new components of the 
same phase until the sum of their Hi column densities is equal to 
or greater than the total Hi column density set for the phase. In the 
latter case, the Hi column density of the last component is reduced 
to a value which makes the sum of Hi column densities match the 
total Hi column density. In this process, if the A'hi of the last com- 
ponent falls below the minimum allowed value (1 x 10" cm"-), all 
the components are discarded and the entire process is repeated. 
Similarly, if the brightness temperature along the line of sight ex- 
ceeds 500 K, then all components are discarded and the process is 
repeated. Note that the highest brightness temperature observed in 
the Galaxy is ^ 200 K, significantly lower than the above threshold. 

The above procedure is carried out for all three phases. Once 
all the components in all three phases have been chosen, the order 
of the components along the sightline is randomized. The known Hi 
column density and spin temperature of each component are then 
used to compute the two observables for the sightline, the total op- 
tical depth Ttot and the brightness temperature Tb; these are then 
used in equations (|4j and {SJ to estimate the Hi column density 
from the optically-thin (^hlot) and isothermal (A'm.iso) methods. 
Of course, the true Hi column density along the sightline, A'hi, is 
already known. For each bin in total "true" Hi column density (bin 
width = 0.235 in log[A'Hi]) and total Hi2Icm optical depth (bin 
width = 0.2), we carried out 2001 runs for each of the six distri- 
butions between the different temperature phases, computing the 
statistics of the ratios (A'hi/A'hi.iso) and (A'hi/A'hi.ot) in every case. 

Finally, we also separately ran the simulations with the tem- 



© 0000 RAS, MNRAS 000, 000-000 



4 Chengalur et al. 



perature ranges 20 < < 200 K for the CNM, 200 K < T, < 
5000 K for the UNM, and 5000 K < T., < 8000 K for the WNM, i.e. 
assuming that the Hi 21cm transition is thermalized (with Ts « Tji) 
in all three phases. This too was done for the above six distribu- 
tions of gas between the different phases. No significant difference 
was found between the results of the simulations using Ts a Tj^ and 
those using the Liszt (^2001) relation between Ts and T^; we will 
hence restrict the discussion to the latter in the following sections. 



4 RESULTS AND DISCUSSION 

In practice, qualitatively similar results were obtained for the dif- 
ferent distributions of gas between the different temperature phases. 
We hence show results only for the last case, where the fraction of 
gas in the CNM, WNM and UNM phases is allowed to vary ran- 
domly from run to run. This is likely to provide the most robust 
estimate of the error on the results, unbiased by a specific choice of 
a gas distribution. 

The critical quantity in determining the accuracy of the 
optically-thin and isothermal estimates of the Hi column density 
is the total optical depth along the sightline Ttot. The accuracy of 
either estimate is worse at high optical depths than for r < 1, be- 
cause it is easier to hide Hi behind foreground gas that is optically 
thick. Figs.[nA] and [B] show, respectively, the ratio of the true 
Hi column density to the optically-thin and isothermal estimates of 
Whi, for the case where the fractions of gas in the three temperature 
phases are allowed to vary randomly from run to run. As expected, 
the optically-thin estimate always yields a lower limit to the true Hi 
column density, i.e. the ratio is always > 1. Even for Ttot ~ L the 
median Hi column density derived from the optically-thin estimate 
underestimates the true Hi column density by a factor of w 1.6; the 
under-estimate is by more than an order of magnitude for optical 
depths of a; 10. The spread on the optically- thin estimate is also 
very large, with the 68.2% confidence level on the ratio A'hi/A'hi.ot 
reaching » 40 for Ttot » 10. 

Conversely, the median isothemal estimate of the Hi column 
density Wkliso tracks the true Hi column density to better than 10% 
even for optical depths of » 5. Even for Ttot ~ 10, the median value 
of Nui/Nnnso is within a factor of 1.5 of unity. The 68.2% confi- 
dence intervals on the ratio extend from » 0.4 to « 3, even at high 
optical depths. Clearly, even for Ttot ~ 10, the isothermal estimate 
of the Hi column density appears accurate to within a factor of » 3. 

Figs.UtA] and [B] show the same ratios plotted in Fig.[T] but 
this time as a function of the true Hi column density. Note that the 
largest total optical depth on sightlines included in this figure is 
» 1000. The left panel [A] shows that, while the optically-thin es- 
timate provides a lower limit to the true A'hi, this is not particularly 
useful at large values of A'hi (<; 10^^ cm"^), for which the true A'hi 
is more than an order of magnitude larger than A'hi.ot- On the other 
hand, the median value of the thin screen estimate A'hi.iso tracks the 
true A'hi to within a factor of 2 over the entire Hi column density 
range (A'hi < 5 X 10^^ cm^^). Further, even the 68.2% confidence 
level intervals of the ratio lie within a factor of a few from unity. 
Clearly, the scatter in the estimated A'hi.iso is quite modest even for 
very high Hi column densities. 

The simulation results thus indicate that the isothermal esti- 
mate of equation ([Sjl provides a fairly good measure of the true 
Hi column density, certainly good to within a factor of a few 
within 68.2% confidence intervals, even for high optical depths, 
Ttot ~ 1000, and high Hi column densities, A'hi ~ 5 x 10^^ cm"-. 
(Note that optical depths of x 1000 are significantly higher than 



would be expected along Galactic sightlines.) This estimator should 
thus allow one to directly probe the existence of very high Hi col- 
umn densities (A'hi <: 10^^ cm"^) in the Galaxy, via accurate mea- 
surements of the Hi 21cm brightness temperature and optical depth 
from, respectively, high velocity resolution Hi 21cm emission and 
absorption spectroscopy. The best fit 2nd-order polynomials to the 
ratio of the true Hi column density to the isothermal estimate of the 
Hi column density are given in the captions to Figs.[T]and|2] Note 
that these are valid for Ttot ^ 1 and A'hi <: 10^' cm"^, respectively. 
Of course, at lower optical depths, Ttot < L the optically-thin esti- 
mate provides an acceptable measure of the true Hi column density, 
within a factor of » 1.5 of the true value. 

It should be emphasized that the isothermal approach as- 
sumes implicitly that the Hi 21cm brightness temperature and op- 
tical depth are measured along the same sightline. Of course, this 
is not the case in reality, with the optical depth along a sightline 
usually measured towards compact background radio sources and 
the brightness temperature inferred by i nterpolating between mea - 
surements at neighbouring positions (e.g.'H eiles & Trolandl2003"ah . 
Small-scale structure in the Hi could imply incorrect brightness 
temperature estimates along the sightline and hence larger errors in 
the estimate of the Hi column density. We do not anticipate that this 
will be a severe problem for most sightlines, especially given the 
possibility of obtaining Hi 21cm emission spectra with telescopes 
of very different angular resolutions (e.g. from the LAB survey, the 
Green Bank Telescope, the Arecibo telescope, etc). 

The approach also assumes that it is possible to accurately 
measure very high Hi 21cm optical depths, Ttot >> 10. In reality, 
direct measurements of such opacities would require spectral dy- 
namic ranges » lO"* per 1 km s"' channel, that will not be achieved 
even with future facilities like the Square Kilometer Array. How- 
ever, it should be possible (except for cases with severe line blend- 
ing) to determine the opacity at line peak via a Gaussian-fitting pro- 
cedure to spectra of a high velocity resolution, as the wings of the 
features (which have lower optical depths) would also have been 
measured at very high signal-to-noise ratios. Note that the high- 
est Hi 21cm opacities would arise for cold gas with high Hi col- 
umn density, which is likely to be thermalized and, hence, to have a 
Gaussian line profile. We have verified from simulations that opaci- 
ties of > 50 can be easily measured via such a Gaussian-fitting pro- 
cedure, from spectra with velocity resolutions of a 0.5 kms"' and 
optical depth RMS noise values of a; 0.001 per 0.5 km s"' channel 
(w hich have already been a chieved with today's interferometers; 
e.g. lBraun & Kanekajl 20051) . 

The present simulations have been carried out so as to match 
observations with high velocity resolution, < 1 kms"'. This is to 
ensure that there is no loss of information due to smoothing of nar- 
row absorption from cold Hi clouds, even at Ti^ a 20 K. The quoted 
errors above on the inferred Hi column density are per 1 kms"' 
spectral channel. However, as noted earlier, the velocity width of 
absorption/emission from individual gas "clouds" with T^ s> 20 K 
will always be significantly larger than 1 km s"' . As a result, the er- 
ror introduced when estimating the total Hi column density along a 
sightline by integrating the per channel isothermal estimate across 
the line profile will be correlated between groups of neighbouring 
1 kms"' channels. A conservative estimate of the net systematic 
error on the total Hi column density along the sightline can be ob- 
tained by assuming that the errors are perfectly correlated over the 
line profile, implying that the error on the total Hi column density is 
the same as that one on the Hi column density per 1 km s"' velocity 
channel, as determined from the present simulation. 

We also note that the assumption of Hi column densities of 
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Figure 1. The ratio of the true Hi column density Nm to the Hi column density obtained from [A] the optically-thin estimate Wrlot (left panel) and [B] the 
isothermal estimate Whi.iso (right panel), plotted against total optical depth Ttot (with a bin width of 0.2). In both panels, the filled circles mark the median 
value of the ratio, while the error bars encompass the 68.2% confidence intervals. For each bin, the median and the co nfiden ce interval were computed over a 
total of 2001 runs. The simulation is for a random distribution of gas between the three temperature phases, using the lLis i3(2001) relation between spin and 
kinetic temperatures. The solid line in the right panel shows the best 2nd-order polynomial fit to the relation between [A'hi/^hi.iso] and riot, valid for Ttoi > 1: 
[A'hi/A'hi.iso] = 0.904 + 0.074rto, - 0.0026r2^,. 
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Figure 2. The ratio of the true Hi column density Whi to [A] the optically-thin estimate Whi.ot (left panel) and [B] the isothermal estimate A'hi.iso (right panel), 
plotted against log[A'Hi] (with a bin width of 0.235 in log[A'Hi]). In both panels, the filled circles mark the median value, while the eiTor bars give the 68.2% 
confidence intervals. For each bin, the median a nd the confid ence interval were computed from 2001 rans. The simulation is for a random distribution of gas 
between the three temperature phases, using the lLiszll l EoOll) relation between spin and kinetic temperatures. The solid line in the right panel shows the best 
2nd-order polynomial fit to the relation between [Whi/Whi.iso] and loglA^Hi], valid for A'hi > 10^' cm"^: [A'hi/A'hijso] = -112.88+ 10. 144 X log [Wni/cm^^]- 
0.2248 xlog[A'Hi/cm"2]-. 



X \Cp - lO^'' cm"^ per 1 kms"' velocity channel is quite extreme. 
Extremely high Hi column densities, > 10^^ cm"^, are more likely 
to be distributed over a significantly larger velocity range, with far 
smaller columns arising per 1 kms"' channel. Allowing for such 
large column densities per 1 km s"' channel results in both shifting 
the median value of the ratio A'hi/A'hi.iso slightly away from unity 
as well as significantly increasing the eiTor bars in the isothermal 
estimate. For example, if we limit to A'hi = 10"^ cm"^ per 1 kms"' 
velocity channel, the median A'hi/A'hi.iso is < 1-2 for r,o, < 10, 
while the 68.2% confidence level intervals on the ratio extend from 
~ 0.5 to ~ 2. Our choice of lO^'' cm"^ as the limiting Hi column 
density per 1 kms"' channel is thus a very conservative one. 

Finally, other estima tors of the Hi column den s ity have been 
used in the literature (e.g. iLockman & Savag3ll995l : IWakker et al.l 



l201ll) . albeit without a detailed characterization of t heir a ccuracy 
in the high A'hi regime. For example, .Wakker et al.l (1201 ll) use an 
estimator that appears similar to the isothermal estimate, with 



1.823 X 10'** X 



/X.log 



Ts - Tb 



dV. 



(10) 



and assume Ts = 135 K for the half of the spectrum with the high- 
est Tb values and T, = 5000 K for the half of the spectrum with 
the lowest Tb values. While this has the advantage that only the 
Hi 2 1 cm emission spectra are needed to determine the Hi 2 1 cm col- 
umn density (as is also true for the optically-thin estimate), it is not 
very useful at high optical depths, where Tj w Tb. It is in the high 
opacity regime that the isothermal estimate is clearly superior to 
the other methods. 
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5 SUMMARY 

We study in detail the "isotiiermal" estimator of the Hi column 
density, earlier proposed by Dickev & Bensoi^ ( Il982h . based on 
measurements of the total Hi 21cm brightness temperature Tb and 
the total Hi 21cm optical depth Ttot from high velocity resolution 
spectroscopy. The "isothermal" estimate is given by A'hi.iso = 
1.823 X 10"* / [Tto, X Tb] / [1 - e""'"'] dV. We have carried out a 
Monte Carlo simulation of realistic sightlines, including gas in dif- 
ferent phases and random locations along the path, to determine the 
accuracy of the isothermal estimate of the Hi column density. We 
find that this approach yields accurate estimates of the Hi column 
density, and that the results do not strongly depend (statistically) on 
the spatial distribution of gas along the sightline, or its distribution 
in different temperature phases. In general, the median isothermal 
estimate of the Hi column density is within a factor of 2 of the true 
Hi column density, while the 68.2% confidence intervals of the ratio 
of A'hi/^hi.iso are within a factor of r; 3 of unity, out to extremely 
high Hi column densities, ^ 5 x 10^' cm"^ per 1 kms"' channel, 
and high total optical depths, < 1000. The 68.2% confidence inter- 
vals are conservative, as they allow for the possibility of extremely 
high Hi column densities, > 10"^ cm"^, in a narrow velocity range 
(1 kms"'), which is unlikely to arise in reality. We conclude that 
the isothermal estimator allows one to accurately measure the Hi 
column density (to within a factor of a few even at the highest Hi 
column densities) and to thus directly probe the existence of ex- 
tremely high Hi column density gas in the Galaxy. 
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